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ABSTRACT 

(D ■ 
pL| I Voids are the most prominent feature of the large-scale structure of the uni- 

^ I verse. Still, they have been generally ignored in quantitative analysis of it, es- 

I sentially due to the lack of an objective tool to identify the voids and to quantify 

^ I them. To overcome this, we present here the VOID finder algorithm, a novel 

lO ! tool for objectively quantifying voids in the galaxy distribution. The algorithm 

I I first classifies galaxies as either wall galaxies or field galaxies. Then, it identifies 

CN I voids in the wall galaxy distribution. Voids are defined as continuous volumes 

j^^ I that do not contain any wall galaxies. The voids must be thicker than an ad- 

^ I justable limit, which is refined in successive iterations. In this way we identify 

'-^ I the same regions that would be recognized as voids by the eye. Small breaches in 

' I the walls are ignored, avoiding artificial connections between neighboring voids. 

H I We test the algorithm using Voronoi tessellations. By appropriate scahng of 

^ I the parameters with the selection-function we apply it to two redshift surveys, the 

'^ I dense SSRS2 and the full-sky IRAS 1.2 Jy. Both surveys show similar properties: 

IT^ \ ~ 50% of the volume is filled by the voids. The voids have a scale of at least 

;h ' 40 h~^ Mpc, and an average under-density of —0.9. Faint galaxies do not fill the 

■ - - ' voids, but they do populate them more than bright ones. These results suggest 

that both optically and IRAS selected galaxies delineate the same large-scale 
structure. Comparison with the recovered mass distribution further suggests 
that the observed voids in the galaxy distribution correspond well to under- 
dense regions in the mass distribution. This confirms the gravitational origin of 
the voids. 



Subject headings: cosmology: observations — galaxies: clustering — large-scale 
structure of the universe — methods: numerical 
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1. Introduction 



Perhaps one of the most intriguing findings of dense and complete nearby redshift sur- 
veys has been the discovery of large voids on scales of ~ 50 h~^ Mpc, and that such large 
voids appear to be a common feature of the galaxy distribution. Early redshift surveys 
like the Coma/Al367 survey (Gregory & Thompson 1978) and the Hercules/ A2 199 survey 
(Chincarini et al. 1981) gave the first indications for the existence of voids, each revealing a 
void with a diameter of ~ 20 h~^ Mpc. Surprising as these findings may have been, it was 
not before the discovery of the Bootes void (Kirshner et al. 1981) that the voids caught the 
attention of the astrophysical community (for a review, see Rood 1988). 

The unexpectedly large void found in the Bootes constellation, confirmed to have a 
diameter of ~ 60 /i^^ Mpc (Kirshner et al. 1987), brought up the question whether the 
empty regions we observe are a regular feature of the distribution of galaxies, or rather rare 
exceptions. Wide-angle yet dense surveys, initially two-dimensional and more recently three- 
dimensional, probing relatively large volumes of the nearby universe, established that the 
voids are indeed a common feature of the large-scale structure (LSS) of the universe. The 
publication of the first slice from the CfA redshift survey (de Lapparent, Geller & Huchra 
1986) introduced the picture of a universe where the galaxies are located on the surfaces of 
bubble-like structures, with diameters in the range 25-50 /i~^ Mpc. The extensions of the 
CfA survey (Geller & Huchra 1989), complemented in the south hemisphere by the SSRS 
and its extension, the SSRS2 (da Costa et al. 1988; 1994) have shown that not only large 
voids exist, but more importantly — that they occur frequently (at least judging by eye), 
suggesting a compact network of voids filling the entire volume. 

It has been recognized early on that inhomogeneities on such scales could impose strong 
constraints on theoretical models for the formation of LSS. However, the voids have been 
largely ignored and their incorporation into our theories of LSS has been relatively recent 
(Blumenthal et al. 1992; Dubinski et al. 1993; Piran et al. 1993). The major obstacle con- 
fronting the incorporation of voids into models of LSS has been the difficulty of developing 
proper tools to identify and to quantify them in an objective manner. As such, the de- 
scription of a void-filled universe with a characteristic scale of 25-50 h~^ Mpc relied solely 
on the visual impression of redshift maps. In order to make a more quantitative analysis 
we have developed an algorithm (void finder) for the automatic detection of voids in 
three-dimensional surveys. Unlike other statistical measures (e.g., the VPF), our target is 
to identify the individual voids, in as much the same way as voids are identified by eye. The 
main features of the algorithm are: 

1. It is based on the point-distribution of galaxies, without introducing any smoothing 
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scale which destroys the sharpness of the observed features. 

2. It allows for the existence of some galaxies within the voids, recognizing that voids 
need not be completely empty. 

3. It attempts to avoid the artificial connection between neighboring voids through small 
breaches in the walls, realizing that walls in the galaxy distribution need not be homo- 
geneous as small-scale clustering will always be present. 

After a review of earlier works (§^, we present a detailed description of our algorithm 
(§0)- We use Voronoi tessellations as a test bed for the algorithm in §^. We then apply the 
algorithm to the SSRS2 redshift survey (§^. Finally (§^, we compare the voids found in 
SSRS2 and in IRAS 1.2 Jy, and discuss the results. 



2. Earlier Works 

The various methods for describing the void content of the LSS of the universe can be 
divided into two categories: statistical measures, and algorithms for identifying individual 
voids within a sample. 

The major statistical tool used for describing the voids is the Void Probability Function 
(VPF). It measures the probability Po{V) that a randomly positioned sphere of volume V 
contains no galaxies (White 1979). For a completely uncorrelated (Poissonian) distribution, 
it is: 

PoiV) = exp{-nV) (1) 

where n is the number-density of galaxies, so that any departure from this quantity represents 
the signature for the presence of clustering. The major drawback of the VPF is that it is 
very sensitive to the details of the galaxy distribution. For instance, adding a few galaxies 
in the under-dense regions may greatly modify the VPF. A less sensitive variant of the VPF 
is the Under-dense Probability Function (UPF), defined as the probability Psp/p{V) that a 
randomly positioned sphere of volume V has a Sp/p under-density (Little & Weinberg 1994). 

The first zero-crossing of the two-point correlation function ^(r) was used by Goldwirth 
et al. (1995) for determining the maximum diameters of voids in case of spherical voids. The 
two-point correlation function is defined as the probability in excess of Poisson distribution 
of finding a galaxy in a volume 6V at a distance r away from a randomly chosen galaxy: 

6P = n6V[l + ^{r)] (2) 
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where n is the mean galaxy number- density. For a cellular like distribution the first-zero 
crossing is a direct measure of the characteristic size of the cells. Using Voronoi tessellations, 
they show that despite the large uncertainty in the determination of ^(r) on large scales 
(> 20 h~^ Mpc), the zero-crossing statistic may be a useful tool in determining the scale of 
typical voids — if the galaxy distribution is void-filled, and if there is a characteristic scale 
for the void distribution. Both these questions are addressed here. Examining the SSRS2 
sample, it was found that -Rzem ~ 38 h~^ Mpc. 

Previous works have used various definitions for voids, and applied different algorithms 
to identify them. Perhaps the first work identifying voids in a quantitative manner is that of 
Pellegrini et al. (1989), who examined ensembles of contiguous cells with densities below a 
given threshold. They use a cubic lattice, and define a local density for each cell in the lattice. 
This local density is based on the analysis of the smoothed density field, constructed from 
the original discrete galaxy distribution. Groups of cells with densities below a specific limit 
constitute the voids. The algorithm considers two cells as contiguous if they are in contact 
either through their faces, edges or vertices. This technique was applied to the original SSRS, 
identifying 4 to 8 voids, depending on the density threshold. Unfortunately, there is only 
partial overlap between the SSRS and the SSRS2. Thus, an overall comparison between the 
results is impossible (where possible, corresponding voids are indicated in Table p. The 
major shortcomings of this algorithm are its use of a smoothed density field, and the lack of 
sense of the shape of the void it recognizes, allowing for practically any void shape. 

Kauffmann & Fairall (1991, hereafter KF91) designed a more elaborate algorithm. They 
too used (empty) cubes, to which adjacent faces are attached. However, in order to avoid 
long finger-like extensions leading from one void into other voids, they impose a constraint 
on the adjacent faces, that each face must have an area of no less than two-thirds that of 
the surface on to which it was to be added. This scheme is restrictive, as it is tailored for 
finding only ellipsoidal-shaped voids. The algorithm was applied to the Southern Redshifts 
Catalog (SRC) and to an all-sky catalog, finding a peak in the spectrum of void diameters 
between 8 and 11 h~^ Mpc. 

Only ~ 10 (out of > 100) of the voids found in KF91 are located within the boundaries 
of the SSRS2. About half of the KF91 voids have an equivalent diameter d < 20 h~^ Mpc, 
which we believe to be statistically insignificant (see § ^.31 ). Other KF91 voids are much 
smaller than their counterparts identified by the void finder: the Sculptor void (void 5 
here) is 40% larger in our analysis, and other voids we find are up to three times larger than 
the corresponding KF91 voids. The limited algorithm and the different statistical analysis 
have resulted in a small void diameter, not describing the true nature of the void distribution. 

A more recent work applying another void search algorithm, is that of Lindner et al. 
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(1995). In this work single spheres, that are devoid of a certain type of galaxies (depending on 
the morphological type and luminosity), are used. The algorithm is applied to an area north 
of the super-galactic (SG) plane, showing that voids defined by bright elliptical galaxies have 
a mean diameter of up to 40 h~^ Mpc, in agreement with the VOID finder results. When 
considering fainter galaxies, the voids are smaller, with the faintest galaxies studied defining 
8 h~^ Mpc voids, suggesting that faint galaxies delineate smaller voids within larger ones, 
which are defined by the bright galaxies. 



3. The Void Finder Algorithm 

The VOID FINDER algorithm was designed with the following conceptual picture in mind: 
The main features of the LSS of the universe are voids surrounded by walls. The walls are 
generally thin, two-dimensional structures characterized by a high density of galaxies. They 
constitute boundaries between under-dense regions, generally ellipsoidal in shape — the voids. 
Although coherent over large scales, the walls — being subject to small-scale clustering — are 
not homogeneous and contain small breaches which we wish to ignore. Galaxies within walls 
are hereafter labeled wall galaxies, while the non-wall galaxies are named field galaxies. The 
voids are not totally empty: there are a few galaxies in them, which we call void galaxies. 
The void galaxies are a sub-population of the field galaxies, indicating which of the field 
galaxies are in the voids. Wall galaxies, on the other hand, cannot — by definition — be in a 
void. 

We define a void as a continuous volume that does not contain any wall galaxies and is 
nowhere thinner than a given diameter. In other words, one can freely move a sphere with the 
minimal diameter all through the void. This definition does not pre-determine the shape 
of the void: it can be a sphere, an ellipsoid, or have a more complex shape, including a 
non-convex one. The definition is targeted at identifying the same regions that would be 
recognized as voids, when interpreting a point distribution by eye. As the voids are defined 
based on the point distribution of galaxies, we do not need to introduce any smoothing 
scale. This is especially important since one of our goals is to measure the volumes of the 
voids. As the regions containing most of the galaxies — walls or filaments — are thinner than 
~ 10 h~^ Mpc, even relatively small smoothing scales will spread the dense regions into the 
under-dense ones, artificially diminishing the voids. 

Our voids may contain galaxies. A stiffer requirement, such that voids should be com- 
pletely empty, is too restrictive as a single galaxy located in the middle of what we would 
like to recognize as a void might prevent its identification. However, for this definition to 
be practical we must be able to identify the field galaxies before we can start locating the 
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voids. 

The algorithm is divided into two steps. First the wall builder identifies the wall 
galaxies and the field galaxies. Then the VOID finder locates the voids in the wall galaxy 
distribution. All together, our method incorporates three parameters, which we define below. 
Two of these parameters (n and /3) are used to determine the field galaxies; the third pa- 
rameter (^) is used during the void search. Specific values for these parameters were chosen 
after a trial-and-error procedure with various simulations (§P, to give results that resemble 
as much as possible eye-estimates of the voids. 



3.1. The Wall Builder Phase 

We define a wall galaxy as a galaxy that has at least n other wall galaxies within a sphere 
of radius i around it. The radius i is hereafter referred to as the wall separation distance. A 
galaxy that does not satisfy this definition is classified as a field galaxy. This is a recursive 
definition which we apply successively until all the galaxies are classified. The minimal value 
that enables filtering of long thin chains of galaxies is n = 3. The choice n = 2 allows thin 
chains extending between dense structures. With n = 3 the minimal structure that can 
be recognized as a wall is a pyramid-like 4-galaxy cluster, with all distances between these 
galaxies smaller than i. The wall separation distance i is chosen in the following manner: 
Let the distance to the n'th closest neighbor of a galaxy be /„. For a given sample this 
quantity has an average value /„ and a standard deviation an- The radius i is defined as: 

i = ln + Pan (3) 

We have chosen here n = 3 and (3 = 1.5. 

Fig. ID demonstrates how the WALL BUILDER works for n = 3. Notice how the galaxy 
string is filtered, while the dense structures (upper and lower parts of each panel) are identi- 
fied and maintained. Additional examples of the wall builder, using Voronoi tessellations, 
are presented in §^. As a side-bonus of this procedure, originally dedicated to filtering the 
field galaxies, we obtain a visual identification of the walls. This is done by drawing all the 
links between wall galaxies satisfying: 

dist{xi,Xj) < i (4) 

These connections are not used in the next step of void analysis, but provide us with another 
visual tool to examine our results (e.g., see Fig. |). 
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3.2. The Void Finder Phase 

The VOID FINDER searches for spheres that are devoid of any wall galaxies. These 
spheres are used as building blocks for the voids. A single void is composed of as many 
(or as few) superimposing spheres as required for covering all of its volume. The algorithm 
is iterative, initially identifying voids containing relatively large empty spheres, than pro- 
ceeding to voids containing smaller spheres. The process can be stopped at any desired void 
resolution, defined as the diameter d of the minimal sphere used as a void building-block. For 
practical reasons the actual process of improving the resolution is done in discrete iterations. 
We denote by di the void resolution used during the i'th iteration. Thus, a void containing 
a sphere of diameter larger than d± would have been identified during the i'th iteration, if 
not earlier. On the other hand, if the largest sphere in a void is smaller than d±, that void 
would not have been identified (yet). 

The motivation for using this gradual process of refining the resolution is the problem 
of keeping apart neighboring voids. If voids were always well-separated by easily recognized 
homogeneous over-dense regions, there would not be any reason to go through this compli- 
cated process. One could simply start from an initial sphere, located somewhere in the void, 
and work his way from that sphere until he is stopped by a wall. Doing this in all directions 
and from all spheres would surely be enough to encompass the whole void. However, the 
situation is more complicated. The walls are not homogeneous, often missing a few bricks 
now and again. Our eyes are capable of ignoring these breaches in the walls, still identifying 
the voids as individual objects. For an algorithm aimed at objectively identifying the voids, 
this is a major problem: how to keep apart neighboring voids, when the boundary is not 
well marked? A simple application of the above process would result in one all-connected 
void, since distinct voids would merge with finger-shaped connections passing through the 
gaps in the walls, in this way connecting those that should have been kept separate. 

The size of the spheres per se is not a good criterion for differentiating between a 
sphere which is a part of the void proper, and a sphere that is passing through a wall 
making a false connection. To illustrate this we will use the somewhat artificial example of 
triangular voids (see Fig. |], left triangle). To cover a large fraction of this void, one sphere is 
perhaps sufficient. But to reproduce faithfully the original void, several spheres with various 
diameters are required — some of the additional spheres being quite small, compared to the 
initial one. Now, imagine a situation where the two voids are not ideally separated — the two 
triangles somewhat overlap (Fig. 0, right-hand side). We find that spheres of the same size 
cover the volume of the voids near the remote edges, but also might connect the two voids 
to one. 

By using a gradually refined void resolution, we overcome this problem in most situa- 



tions. At the risk of being tautological, we will note that the difference between the void 
proper and a false connection is that the connections are characterized by them being thinner 
than the rest of the void and leading into another void. These two properties are critical for 
our purpose. The first property allows for the initial identification of the center of the void, 
before finding out about the connections. Using a relatively crude resolution, we will first 
find the larger spheres that actually make up the void. Now, when the resolution becomes 
sufficiently small to allow spheres to cross the walls, we will discover that the other side of 
the wall is already occupied by another void (identified earlier). 

Hence, the spheres for covering a void are picked up in two stages: the identification 
stage, followed by consecutive enhancements. We will now describe these stages in detail.^ 

The identification stage identifies the central parts of the void. Usually, these spheres 
cover only about half of the actual volume. We focus (at this stage) on identifying a certain 
void as a separate entity, rather than trying to capture all of its volume. The central parts 
of a void are covered using spheres with diameters in the range ^d^ax < d < d^ax, with 
c^max denoting the diameter of the void's largest sphere. The parameter ^ is the thinness 
parameter. Once a group of such intersecting spheres has been dubbed a void, it will not be 
merged with any other group. The parameter ^ controls the flexibility allowed at this stage. 
Setting ^ = 1 would leave us with only the largest sphere in the void, while lowering ^ allows 
the addition of more spheres. If the void is composed of more than one sphere (as is usually 
the case), then each sphere must intersect at least another one with a circle wider than the 
minimal diameter C,dma,x- We have chosen ^ = 0.85, which allows for enough flexibility — 
without accepting counter-intuitive void shapes. A lower ^ reduces the total number of the 
voids, with a slow increase in their total volume. 

After the central part of a void has been identified, we consecutively enhance its volume, 
in order to cover as much of the void volume as possible using the current void resolution. 
These additional spheres need not adhere to the ^ thinness limitation: we scan the immediate 
surroundings of each void, and if empty spheres are found then they are added to the void. 
We scan for enhancing spheres of a certain diameter only after scanning for new voids with 
that diameter. In this way we do not falsely break apart individual voids, and we do not 
prevent the identification of truly new voids. By applying consecutive enhancements to the 
voids, we have elaborated our definition for a void, by allowing the voids to gradually become 
thinner. Therefore the precise void definition is a continuous volume that does not contain 



^This later addition of spheres to the central part of the voids is the main improvement here over the 
VOID FINDER algorithm that was presented earlier (El- Ad, Piran & da Costa 1996a). The older version, 
initially used for analyzing the SSRS2, did not include this feature. All the results presented here use the 
improved version of the algorithm. 
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any wall galaxies, and is thicker than an adjustable limit. If this definition seems complex, it is 
a tribute to the difficulty of teaching a machine how to grasp a composite three-dimensional 
object. 

Using once more the triangular voids example (Fig. ^, right hand side), notice that we 
do not falsely identify a new void near one of the edges of the triangle, since the initial 
sphere gradually takes-up that volume, leaving room only for smaller spheres which are well 
below the current void resolution. Thus we identify only one void per triangle, as we should. 
Further still, we do not connect the two voids, as they are initially identified as two separate 
entities — not to be merged later. 

To conclude, we present a step-by-step listing of the various stages: 

1. Choose (ii. 

2. Find all empty spheres with d > C,di. 

(a) Construct a grid with a cell size {$,d±/^/3Y. 

(b) Choose the centers {cj} of the empty grid cells that are not contained within 
previously found voids. 

(c) Starting from each Cj, construct a continuous path along which the radius of the 
maximal empty sphere centered on the path continuously increases. The path 
ends at a point c™*^^ where the size of the maximal sphere is a local maximum. 

(d) Keep all maximal spheres having d > C,di. The redundancy in the grid size insures 
that every sphere with d > C,d± contains at least one Cj, so the identification of all 
such spheres is guaranteed. 

3. Arrange the spheres into groups (i.e., new voids), satisfying: 

(a) Each group contains at least one sphere with d > d±, the largest of which is 
denoted (imax- 

(b) A group may include additional spheres; these intersect at least one other sphere 
in the group with a circle whose diameter is larger than or equal to ^cimax- 

4. Enhance the older voids: 

(a) Scan the outskirts of all previously identified voids for empty spheres. 

(b) Add d > di spheres that intersect a void to it. 

5. Decrease the void resolution d±, and iterate the process. 
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Fig. ^ demonstrates how the void finder works. In all panels we show the same 
slice through a three-dimensional Voronoi Tessellation (see §^ for details), and we follow the 
development of the void image as we refine the void resolution. Initially, we locate the voids 
containing the largest empty spheres. In the following iterations we locate the smaller voids, 
and — when appropriate — enlarge the volumes of the older ones. 

One must take into account the boundaries of the distribution. We treat the limits 
of the distribution as rigid boundaries. This may cause some distortion in the voids found 
close to the boundaries (causing them to be smaller than the rest). Still, this is the least 
speculative option as we do not make any assumptions about unsurveyed regions. 



3.3. Statistical Significance 

To assess the statistical significance of the voids we compare the voids found in observed 
data with voids found in equivalent random distributions. The random distributions mimic 
the sample's geometry and density, and are analyzed by the algorithm in exactly the same 
manner. Averaging over the random catalogs we calculate iVpoisson(c?), the expected number 
of voids in a Poisson distribution as a function of the void resolution d. We compare this 
with the observed number, iVsurvey('^)- We define the confidence level as: 

P{d) = 1 - -T-r TTT (5) 

^Vsurveyl"j 

The closer p{d) is to unity, the less likely the void could appear in a random distribution. 
We consider voids with p > 0.95 as statistically significant. At a certain void resolution t/gtop, 
A^Poisson('^) exceeds iVsurvcy('^)! and we terminate the void search. 

The statistical significance defined in this way is attributed according to the order in 
which the voids were identified, not according to their sizes. This is so because the parameter 
determining how early a void is detected is the diameter of the largest sphere it contains, 
and not its total volume. Of course the total volume of a void and the diameter of the 
largest sphere contained in it are correlated. However, a void composed of a single large 
sphere will be detected earlier — and hence considered more significant — than a larger void 
composed of several spheres all having smaller radii. This is in agreement with the fact that 
clear spherical voids are more prominent when a galaxy distribution is inspected by eye. In 
retrospect this choice is also justified by the theoretical expectation that voids become more 
spherical with cosmological time (Blumenthal et al. 1992). 

Voids in the actual surveys are systematically larger than the voids in the random 
distributions. When we reach the void resolution dgtop (where the number of voids in the 
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random distributions exceeds the void number in the actual survey), the volume taken up 
by the random voids is still much smaller than the volume of the true voids. By ignoring 
this fact we actually under estimate the significance of the voids. 



3.4. Redshift Surveys 

The average galaxy number- density decreases with depth in a magnitude- limited redshift 
survey. If not corrected, this selection effect will interfere with the algorithm in the deeper 
regions of the sample: field galaxies will occur more frequently, and the derived size of 
the voids will be larger. Consequently, systematically larger voids will be found at greater 
distances. 

To avoid these effects, one should use a volume-limited sample, in which the galaxy 
number-density is constant and independent of the distance. A volume-limited sample with 
M < Mo has a depth: 



where rub is the survey's magnitude limit and Co is the luminosity that corresponds to Mo. 
However, current volume-limited samples are too small to study the LSS. For example, the 
SSRS2 sample contains 3162 galaxies with r < 130 h~^ Mpc. A volume-limited sample based 
on this distribution with such a depth {Mo = —20.1) retains only 528 galaxies. To overcome 
this, we use a semi-volume-limited sample: volume-limited up to some medium radius Tq, 
and magnitude-limited beyond. We choose the depth Tq by maximizing the number of bright 
galaxies N{M < Mo): 

N{M < Mo) = y r/ ■ r]T{l - a, Xo) (7) 

where rj is the galaxy number-density. The incomplete F-function arises from the integration 
of the appropriate Schechter function, with x = Cj L^,. 

In Fig. ^we plot the normalized N{Mo) for the SSRS2 parameters: a = 1.2, M^ = —19.5 
and nib = 15.5. In this example N peaks close to M^ at M = —19.1, which corresponds 
to Xo = 0.69 and a depth Vo = 83 h~^ Mpc. If we increase the depth of the volume-limited 
sample, the number of galaxies decreases. 

No corrections are needed in the volume-limited region. We determine the values for i 
and d± in this region. Beyond To we define 0(r), a selection-function based on the Schechter 
luminosity- function: 

T{l-a,XM) . . 

nr) = T^^ r (8) 

F(l -a,XMj 
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where xm = lO~^-'^('^~'^*\ The selection-function 0(r) is the observed fraction of galaxies at 
the distance r, relative to Tq. A plot of the selection-function for the semi-volume-liniited 
sample for SSRS2 is shown in Fig. ^ Using the selection-function, we modify both phases 
of the algorithm. In the WALL BUILDER phase, we consider larger spheres when counting 
neighbors at r > Vo- The radius of the counting-spheres is modified: 

[ £/(f)^'^(r) otherwise 

A similar correction is applied to the void finder phase. A void of a given size found 
in a low density environment is less significant than a void of the same size found in a 
high density environment. In order that all the voids found in a given iteration are equally 
significant, we adjust the algorithm so that at a given iteration relatively larger voids are 
accepted, if located at r > Vo'. 

< = l^\.i/3. ^ '"l^"- (10) 

[ di/(p ' [r) otherwise 

As mentioned earlier ( §3.3| ), we use random distributions that mimic the true sample's 
geometry and density, in order to assess the statistical significance of the voids. This scheme 
needs to be fitted, to correspond to the decrease in the galaxy density with r. We continue 
to use the same definition for the confidence level p (eq. 0). However, d is now a function 
of r, scaled in the same manner given by equation (|TD|) (see Fig. ^. 



4. Voronoi Distributions 

As a test-bed for the void finder algorithm, we use Voronoi distributions: A dis- 
tribution of galaxies that is based on a Voronoi tessellation (Voronoi 1908). A Voronoi 
tessellation is a tiling of space into convex polyhedral cells, generated by a distribution of 
seeds. To generate a galaxy distribution in which the galaxies are located on the walls of 
the Voronoi cells, we have used an algorithm developed by Van de Weygaert & Icke (1989). 
The resultant galaxy distribution has the desired characteristic of large empty regions (i.e., 
voids), which we identify by the void finder algorithm. 

A Voronoi tessellation is constructed from a given set of seeds {xi}. Based on the 
locations of these seeds, we divide the volume into cells. The Voronoi cell IIj of seed i is 
defined by the following set of points x: 

IIj = { X I dist(x, Xi) < dist(x, Xj) for all j ^ i} (11) 
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In other words, IIj is the set of points that is nearer to the seed Xi than to any other seed. 
In Fig. ^ we show two-dimensional cuts through a Voronoi tessellation. 

We assign a finite width to the walls and position the galaxies on the boundaries between 
the Voronoi cells, with a Gaussian displacement in the distance from the exact cell boundary. 
These galaxies are designated to be wall galaxies, but do not necessarily end-up in this way: 
in regions where the galaxy distribution is sparse, the wall builder may identify these 
galaxies as field galaxies. Further still, if the galaxy density is low, there may be boundaries 
between Voronoi cells with no galaxies at all. Additional random galaxies correspond to field 
galaxies. Most of these randomly placed galaxies end-up as field galaxies. A few, located in 
dense areas, will be identified by the wall builder as wall galaxies. 

We will call a galaxy distribution constructed in this way a Voronoi distribution. The 
location and number of the Voronoi cells (the would-be voids), the spread of the wall galaxies 
and the fraction of random galaxies are all known. Therefore, we can use this distribution 
as a test bed for our algorithm. 

We have constructed a Voronoi distribution using 3000 galaxies and 10 seeds, which 
includes 300 (10%) random galaxies. These numbers correspond roughly to what is available 
in present redshift surveys. The original Voronoi tessellation (Fig. ^ is compared to the 
VOID FINDER reconstruction (Fig. ^. All Voronoi cells are reproduced except the very small 
cells near the boundaries, that were cut by the box limit. The reconstructed voids follow 
closely the original Voronoi cells, withstanding the noise introduced by the random galaxies. 
The walls highlighted using the wall builder are located along the boundaries between 
the Voronoi cells. 

It is worthwhile to examine where the algorithm fails. A single void was broken into two 
(Fig. ^, lower-left corner of upper-left box) due to some random galaxies that extended the 
walls into this cell. Also note the instances where two Voronoi seeds are situated relatively 
close to each other (e.g., in the lower-left box). In such a situation the resultant wall is very 
sparse. This is due to the way in which the Voronoi code generates the galaxy distribution. 

We have also created mock surveys, based on Voronoi distributions. Galaxies in the 
Voronoi distribution were assigned magnitudes according to a Schechter function. Then, a 
magnitude-limited sample was chosen. To this Voronoi-based mock survey we have applied 
our usual procedure: analyzing a semi-volume-limited sample, applying corrections beyond 
To- Fig. ^ shows a reconstruction based on the original data, along with a reconstruction 
based on the corresponding mock survey. The fit between the Voronoi cells and the recovered 
voids is still good, showing the adequacy of our method in analyzing actual surveys. 

All together, the Voronoi tessellations that we examined show that the void finder 
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indeed generates a faithful reproduction of the Voronoi cells. Further still, in cases where 
the reproduction merges adjacent Voronoi cells into one void, we see this as the adequate 
outcome of a missing wall. If we would have examined such a galaxy distribution by eye, 
with no prior knowledge about the locations of the Voronoi cells, we too would most likely 
consider that volume — originally occupied by two Voronoi cells — as one void. A level of 
~ 10% random galaxies is tolerated, with no significant distortion in the void reproduction, 
and Voronoi-based mock surveys are also reproduced faithfully. 

Boundary distortions are the cause for most of the cases in which the void finder 
departs from the original tessellation. This is also evident when considering the void volumes: 
The volume occupied by the voids is ~ 15% larger, if we consider only an inner cube and not 
the full tessellation. Voids near the boundaries are typically smaller (if at all recognized), as 
the Voronoi cells are bisected by the cube's boundary. 



5. Voids in the SSRS2 

The SSRS2 survey (da Costa et al. 1994) consists of ~ 3600 galaxies with mi, < 15.5 
in the region —40° < 5 < — 2?5 and h < —40°, covering 1.13 sr. We have considered a 
semi-volume-limited sample, in this case consisting of galaxies brighter than Mg < —19, 
corresponding to a depth r^ = 79.5 h~^ Mpc. The Schechter luminosity function was evalu- 
ated with M* = —19.5 and a = 1.2, as derived for the SSRS2. Our final semi-volume-limited 
sample consists of 1898 galaxies, extending out to rmax = 130 h^^ Mpc where the selection- 
function (j) has dropped to 17% (Fig. H). It should be emphasized that the SSRS2 analysis is 
performed in redshift- space. However, because of the paucity of large clusters and the small 
amplitude of peculiar motions in the volume surveyed by the SSRS2, redshift distortions are 
small (da Costa et al. 1997) and the properties derived here should reflect those of voids in 
real-space. 

The WALL BUILDER analysis of the SSRS2 has classified 1736 galaxies (91.5%) as wall 
galaxies, and 162 (8.5%) as field galaxies. The wall separation distance was I = 7.4 h~^ Mpc. 
In the volume-limited region we have n~^^^ = 6.4 /i~^ Mpc, so i/n~^^^ = 1.16. The wall 
galaxies are grouped in ten structures. The remarkable fact that one structure contains 
most (96%) of the wall galaxies will be discussed elsewhere, as this may have interesting 
statistical implications regarding the connectivity of the wall-filament skeleton. The rest of 
the wall galaxies are found in nine groups, each having 4 to 21 galaxies. 

We have identified eleven significant {p > 0.95) voids within the volume probed by the 
SSRS2. These voids were detected while the void resolution was d > 19.2 h^^ Mpc. In 



-15- 



the following calculations we take into account only these voids, unless otherwise specified. 
Seven additional voids were identified before the void search was terminated at the resolution 
c^stop = 15.1 h~^ Mpc, for which p vanishes. Fig. ^ plots the accumulated number of voids in 
the SSRS2 and in the corresponding random catalogs, as a function of the void resolution 
d. The locations and characteristics of all eighteen voids are given in Table |l|. Column (1) 
lists the confidence level p attributed to each void. The diameters given in column (2) are of 
a sphere with the same volume as the whole void, as is listed in column (3). The center of 
the void given in columns (4)-(6) is defined as its center-of-(no)-mass. The density contrast 
value listed in column (7) was corrected for the average galaxy density at the same distance 
as the center of the void. In column (8) we give the fraction of the total volume of the void 
covered by the single largest sphere contained in it. This value is typically ~ 40% of the total 
volume of the void. Finally, in column (9) we indicate the corresponding voids identified in 
da Costa et al. (1988, hereafter dC88). 

The average size of the voids in the survey as estimated from the equivalent diameters 
is c? = 40 ± 12 /i~^ Mpc. The average under-density within the voids was found to be 
5p/ p ~ —0.9, a quite remarkable result showing how empty voids are of bright galaxies. The 
eleven significant voids comprise 54% of the survey's volume. An additional 5% is covered by 
the seven additional voids, totaling in ~ 60% of the volume being occupied by these voids. 
We estimate that the walls occupy less than 25% of the volume. 



The SSRS2 sample and the resultant VOID finder reconstruction is presented in Fig. |10, 
as six 6?25-wide constant-declination slices. Do not be mislead by the wide right-ascension 
range of this survey, emphasized by these images. In declination the survey spans 37?5 (see 



Fig. |Tl|), limiting the sizes of the voids we find. The slices contain numerous walls and 
filaments. The nearest prominent features, the Fornax cluster and the Eridanus group, can 
easily be recognized at (r ~ 10 h~^ Mpc, a = 3^5) running through most of the declination 
range. One of the outstanding structures in this survey is the Southern Wall (SW). This 
wall (dC88) is the dense structure running across all slices, from {a = 4^, r = 45 h~^ Mpc) 
to (a = 0^, r = 90/i~^Mpc). It is especially prominent in the 6 > —15° and the 6 < 
— 27?5 slices. The Pavo-Indus-Telescopium (PIT) supercluster runs along the line of sight at 
a = 21^5 in the 6 < — 21?3 slices. These two walls are almost parallel, and are connected 
through a third wall running perpendicular to the line of sight at r ~ 90 h~^ Mpc in the right 
ascension range 21^5 to 0^. These walls bound two voids (void 5 and void 12 in Table |l|), 
identified earlier as one void (void 3 in dC88) — here the void is broken into two by a filament 
almost connecting PIT and the SW. A third structure runs parallel to the SW, bounding 
another void (void 1) between them. 

The largest void found in the SSRS2 survey (void 3) has an equivalent diameter d = 
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60.8 h ^ Mpc, making it comparable in volume to the large void found in the Bootes (Kirshner 
et al. 1981). This void is an ellipsoid, whose major axis is perpendicular to the line of sight, 
located at: 80 h'^ Mpc < r < 130 h-^ Mpc; -25° <6 < -2?5; 21^ < a < 23'?75. This void 
might actually be larger, since it is bounded (in three directions) by the limits of the SSRS2. 
A second large void (void 2, with d = 56.2 h~^ Mpc) is also comparable to the Bootes void. 

When preparing the semi-volume-limited sample, we cast off all faint M > Mo galaxies 
in the region r < r^. These galaxies comprise the bulk of the surveyed galaxy population 
that we are forced to ignore (the rest are r > rmax galaxies, where the sample is too sparse). 
Although we cannot use these galaxies during the analysis phases, we can still try and benefit 
from them a posteriori: after the voids are located, we examine the locations of these galaxies. 
In this way we have divided the original galaxy distribution in the volume-limited region into 
two distinct populations, according to some absolute-brightness limit. 

There are 1264 Mo > — 19 galaxies in the volume-limited region of the SSRS2. Almost 
61% of this region is covered by voids — but only 19% of the faint galaxies are found within 
them. Even though the VOID FINDER algorithm uses only the brighter galaxies in this region, 
we find that the faint galaxies do not fill the voids, providing an excellent verification of our 
algorithm. Still, the percentage of faint galaxies within the voids is significantly larger than 
that of the bright galaxies: only ~ 5% of the bright Mo < — 19 galaxies are contained in the 
voids. 



6. Discussion 

We have used the VOID finder algorithm to analyze two redshift surveys: The SSRS2 
and the IRAS 1.2 Jy survey (El- Ad, Piran & da Costa 1996b). These surveys represent 
two extreme cases in the trade-off between density and sky coverage. The SSRS2 is densely 
sampled (m^ < 15.5), but it has narrow angular limits, especially in the declination range. 
As a result, voids are often limited by the survey's boundary, diminishing their scale. The 
IRAS is an almost full-sky survey (87.6% coverage), but it is rather sparse. As a result one 
cannot use a small void resolution, for lack of statistical significance. In addition to the above 
differences, the SSRS2 galaxies are optically selected, as opposed to the IRAS galaxies. 

Withstanding these differences, the results obtained with these surveys are similar. 
First, the surveys agree regarding individual voids in the regions where the surveys overlap. 



Fig. |T2| depicts the redshift-space voids in the SG plane, for the IRAS and for the correspond- 
ing part of the SSRS2. In the region where the SSRS2 sample overlaps the IRAS sample, we 
find three of the eleven significant voids identified in the SSRS2. The corresponding IRAS 
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voids are ~ 33% larger the SSRS2 ones, since they are not bounded by narrow angular limits 
as the SSRS2 voids. 

Additionally, the results agree viz a viz the statistical characteristics of the voids: 

1. Large voids occupy ~ 50% of the volume. 

2. Walls occupy less than ~ 25% of the volume. 

3. A void scale of at least 40 h^^ Mpc, with an average under-density of —0.9. 

4. Faint galaxies do not "fill the voids", but they do populate them more than bright 
ones. 

The void scale derived in both surveys is a lower limit: for the SSRS2, because of the 
narrow boundaries limiting the voids; and for the IRAS, due to the larger d^topi and because 
of the conservative analysis applied regarding the ZOA and the field galaxies. Further 
still, regarding the walls, both surveys poses the remarkable characteristic where almost all 
(~ 95%) of the wall galaxies are contained in a single structure. 

The fact that both the IRAS and the SSRS2 are consistent regarding the void statistics 
as well as the individual voids is not trivial, since the IRAS galaxies represent a special 
galaxy class, which is biased relative to the optical galaxies (Lahav, Rowan-Robinson & 
Lynden-Bell 1988; Lahav, Nemiroff & Piran 1990). The agreement between the surveys 
suggests that a similar void scale exists for both optically and IRAS selected galaxies. 

The IRAS data has been used to derive the smooth density field, and it probes a volume 
comparable to that used to determine the density field of the underlying mass distribution 
from the POTENT reconstruction method (Dekel, Bertschinger & Faber 1990) based on the 
measured galaxy peculiar velocity field. The voids and walls identified by our algorithm 
(El-Ad, Piran & da Costa 1996b) indeed correspond to the under-dense and over-dense 
regions in the IRAS density field (Strauss & Willick 1995) respectively. Comparison with 
the SFI sample (da Costa et al. 1996) also demonstrates that the voids delineated by galaxies 
correspond remarkably well with the under-dense regions in the reconstructed mass density 
field derived from peculiar velocity measurements (but also compare with the Mark III 
map— Dekel 1994; Dekel et al 1997). 

Most of the over-dense regions, walls and filaments, are narrower than 10 h^^ Mpc. The 
smoothing scale used for creating the density fields spreads the originally thin structures over 
wider regions, extending into the under-dense volumes. This has the effect of giving a false 
impression of a rather blurred galaxy distribution, where prominent over-dense structures 
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are separated by small under-dense regions. The true picture is very different: there is a 
sharp contrast between the thin over-dense structures which occupy only the lesser part of 
the volume, and the large voids. The notion of a void filled universe cannot be avoided in 
this picture. 

We have developed and tested a new tool for quantifying the large-scale structure of the 
universe. We focus on the under-dense regions, and for the first time we have individually 
identified and statistically quantified the voids. The VOID finder analysis clearly shows the 
prominence of the voids in the LSS, not hindered by smoothing of the over-dense regions, 
and it reveals the image of a void-filled universe, where large voids are a common feature. 

The consistency in the void image between IRAS and optically selected galaxies suggests 
that galaxies of different types delineate equally well the observed voids. Therefore galaxy 
biasing is an unlikely mechanism for explaining the observed voids in redshift surveys. Com- 
parison with the recovered mass distribution further suggests that the observed voids in the 
galaxy distribution correspond well to under-dense regions in the mass distribution. This 
confirms the gravitational origin of the voids (Piran et al. 1993). 

We would like to thank Luiz da Costa for providing us with the SSRS2 data and for 
numerous helpful discussions. 
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Confidence 


Equivalent 


Total 


Location of Center 


Void 


Largest 






Level 


Diameter 


Volume 


(Equatorial Coordinates) 


Under- 


Sphere's 






P 


[h-iMpc] 


[/i-3 KMpc^] 


r 


a 


5 


dcnsity 


Fraction 


Identification 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


>0.99 


54.3 


84.9 


85.7 


l'^33™ 


-16°45' 


-0.89 


0.25 




2 


>0.99 


56.2 


93.9 


99.7 


3 34 


-28 50 


-0.87 


0.20 


dC88 void 4 


3 


>0.99 


60.8 


119.0 


107.2 


22 25 


-14 46 


-0.93 


0.25 




4 


>0.99 


35.6 


24.0 


66.7 


21 43 


-14 40 


-0.91 


0.37 




5 


>0.99 


34.8 


22.6 


53.0 


23 48 


-24 39 


-0.94 


0.39 


dC88 void 3 (Sculptor) 


6 


>0.99 


32.0 


17.4 


56.5 


3 56 


-20 11 


-0.92 


0.47 




7 


>0.99 


25.5 


8.8 


77.2 


3 17 


-11 40 


-0.91 


0.71 




8 


>0.99 


27.8 


11.4 


83.9 


23 20 


-12 32 


-0.95 


0.59 




9 


0.99 


39.0 


31.1 


114.6 


3 06 


-13 47 


-0.86 


0.54 




10 


0.95 


34.8 


22.4 


104.7 


26 


-9 17 


-0.69 


0.39 




11 


0.95 


42.9 


41.5 


112.8 


21 


-29 43 


-0.88 


0.31 




12 


0.72 


25.0 


8.1 


74.8 


23'^03™ 


-32°35' 


-0.97 


0.37 




13 


0.72 


22.1 


5.8 


31.0 


1 23 


-19 36 


-1.00 


0.55 


dC88 void 1 


14 


0.53 


21.3 


5.2 


87.2 


21 28 


-29 28 


-1.00 


0.62 




15 


0.53 


27.3 


10.7 


116.1 


21 24 


-33 17 


-0.74 


0.86 




16 


0.36 


20.3 


4.4 


36.5 


21 43 


-18 41 


-1.00 


0.52 




17 


0.27 


19.0 


3.7 


32.1 


3 42 


-21 21 


-1.00 


0.62 




18 


0.27 


21.1 


4.9 


85.9 


4 18 


-8 42 


-1.00 


0.50 





Table 1: Locations & properties of the voids in the SSRS2 survey. 
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Fig. 1. — Wall construction using the WALL BUILDER. Panel a: A toy distribution of 16 
galaxies (o). Panel b: After the calculation of i, all galaxy pairs closer than this separation 
are marked. Panel c: Galaxies with less than three neighbors are flagged as field galaxies 
(©). Panel d: The final result (after one additional iteration of neighbors count). The string 
extending between the dense structures has been eliminated. 
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Fig. 2. — Void coverage using spheres: Larger spheres (hght gray) are identified first, and 
smaller ones (darker gray shades) are identified later. Left: A triangular void. Seven spheres 
cover the void, the radius of the smallest spheres is a quarter of that of the largest one. 
Right: Two triangles partially overlapping, as if there is a breach in the wall separating the 
voids. The gap is larger than the smallest sphere used to cover each void. However, we still 
can separate the two voids. 
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(a) Voronoi seeds and cells (b) ...with the galaxies 



(c) ...& the walls 



(d)d = 0.816 




(e) d = 0.680 



(f)d = 0.510 



(g)d = 0.410 



(h) d = 0.356 



Fig. 3. — A demonstration of the way the VOID FINDER covers the voids. All panels depict 
the same slice, cut through a certain Voronoi tessellation. We present the Voronoi seeds (□) 
& cells (panel a), the galaxies (panel b) and the walls (panel c). The remaining panels show 
the voids' image, at various void resolutions d^^. More voids are recognized as we refine d^, 
and the older voids are enlarged. 
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Fig. 4. — The relative number A^/A^max in the SSRS2 sample {rrih = 15.5, a = 1.2 and 
M* = —19.5) as a function of the depth of the volume-limited region. The number of 
galaxies A^ peaks at M = -19.1 (r = 83 h'^ Mpc). 
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Fig. 5. — Nine planar slices showing Voronoi cells and galaxies in cubic Voronoi tessellation, 
generated from 10 seeds with 3000 galaxies (10% random). The cells are depicted using 
various shades of gray, indicating the intersection of the Voronoi cell with the plane at the 
center of the slab. The locations of the Voronoi seeds are marked by '□'. Galaxies associated 
with cell boundaries are marked by 'o', and random galaxies by '•'. 
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Fig. 6. — The void finder reconstruction corresponding to Fig. |^: The voids are indicated 
using various gray shades, where the depicted voids correspond to the intersection of the 
central plane of each slab with the three-dimensional voids. Also shown are the walls (dark 
lines — connections between nearby galaxies which satisfy eq. [Q) and the galaxy breakdown 
(wall galaxies: 'o' or '•'; field galaxies, outside the voids: '©'; void galaxies: '©'). Non-void 
galaxies that look as if they are in the voids are actually not — they only seem so, due to 
projection effects. Voids in the first and last panels are smaller due to boundary effects. 
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Fig. 7. — Left: void finder reconstruction of an ordinary Voronoi distribution. Right: 
reconstruction of a mock survey based on the same data. The intermediate arc marks the 
volume-hmited region of the mock sample. 
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Fig. 8. — The SSRS2 selection-function, with Vo = 79.5 h^^ Mpc and rmax = 130 /i^^ Mpc. 
The smooth hne depicts 0, where as the dashed hne is the actual relative density, evaluated 
at 10 h^^ Mpc intervals. 
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Fig. 9. — The accumulated number of voids as a function of the void resolution d, for the 
SSRS2 and for equivalent random catalogs. The derived confidence level p is also indicated. 
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Fig. 10. — LSS in the SSRS2. The survey spans the dechnation range —40° < 5 < — 2?5, 
which we have spht to six 6?25-wide shces. Our sample extends to rmax = 130 h~^ Mpc, and 
the intermediate arc marks r^ = 79.5 h^^ Mpc {Continued on the following page). 
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Fig. 10. — (Continued from previous page) Voids (grey areas), walls (dark lines) and galaxies 
(wall galaxies: 'o'; void galaxies: '0'; field galaxies: '©'). 
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Fig. 11. — A constant right-ascension slice of the SSRS2, demonstrating the narrow (37?5) 
dechnation range of this survey. Most of the voids in this image are hmited by the survey's 
boundary. 
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Fig. 12. — Voids in the SG plane: here we compare the redshift-space voids of the SSRS2 
{left) and the IRAS (right). The denser samphng of the SSRS2 is evident. Similar voids are 
found in the overlapping regions of the two surveys. 



